First-principles calculations of spin spirals in Ni2MnGa and Ni2MnAl 
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We report here non-coUinear magnetic configurations in the Heusler alloys Ni2MnGa and Ni2MnAl 
which are interesting in the context of the magnetic shape memory effect. The total energies for 
different spin spirals are calculated and the ground state magnetic structures are identified. The 
calculated dispersion curves are used to estimate the Curie temperature which is found to be in good 
agreement with experiments. In addition, the variation of the magnetic moment as a function of the 
spiral structure is studied. Most of the variation is associated with Ni, and symmetry constraints 
relevant for the magnetization are identified. Based on the calculated results, the efi'ect of the 
constituent atoms in determining the Curie temperature is discussed. 



I. INTRODUCTION 

Materials showing strong coupling between the mag- 
netic and structural properties are interesting from the 
technological point of view. Tb-Dy-Fe alloys (Terfenol- 
D, already in commercial use ) exhibit magnetic field 
induced strains q£ ^ 0.1 % based on the magnetostric- 
tion phenomenorc. On the other hand, Ni-Mn-Ga alloys 
close to the Ni2MnGa stoichiometry-sliow strains up to 
10 % with moderate magnetic fieldad'Ela. The mechanism 
of this phenomenon, the magnetic shape memory (MSM) 
effect, is based on the magnetic-field-induced movement 
of structural domains (twin, variants) and is different from 
ordinary magnetostrictionS. The basic magnetic proper- 
ties related to the MSM effect include the saturation mag- 
netic moment and the magnetic amsotropy which have 
been studied earlier for Ni2MnGaQl3. Here, we probe 
deeper into the magnetic properties of Ni2MnGa and 
an other MSM candidate, Ni2MnAl, by studying non- 
collinear magnetic configurations which also enables one 
to consider finite temperature effects in a natural way. 

Although one ingredient in the MSM effect is a struc- 
tural transformation (martensitic transformation) from 
a cubic structure to a lower-symmetry structure upon 
cooling, we concentrate here only in the high tempera- 
ture phase. In this phase Ni2MnGa has the cubic L2i 
structure (see Fig. |l|) as-|Shown by x-ray and neutron 
diffraction measurementsoB. The magnetic order is fer- 
romagnetic, and most of the magnetic moment originates 
from MnErB. In the stoichiametric compound the Curie 
temperature is about 370 KlJ and decreases when increas- 
ing the Ni contentliy. On the other hand, Ni2MnAl is 
less studied and its structure and magnetic configura- 
tion do not seem to be perfectly understood. On the 
structural sidfty-|bath| J*2i and disordered B2 structures 
are reportedMO't^loEj depending on the thermal treat- 
ment. The magnetic configuration is found to be ferro- 
magnetic with Curie temperatures between 300 K and 
400 K in Ref. |l^ and antiferromagnetic or spiral in Refs. 
|ri|,|l2| The magnetic momen^-nomes mainly from Mn 
atoms also in this compoundOlla. It seems that the 



ground state magnetic configuration depends on the un- 
derlying crystal structure. Here we address the possibil- 
ity of non-coUinear magnetic configurations in the L2i 
structure. 
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FIG. 1: Cubic cell of the L2i structure, where X is Al or Ga. 
The cubic cell contains four primitive cells. 

Although the origiuaj. formulation of the local-spin- 
density approximationliZI of density-functional theory al- 
lowed non-coUinear magnetic order, first-principles cal- 
culations for this aspect have started only recently (for a 
review see Ref. |l^). One application has been the-studji 
of non-coUinear ground states for |e«jajaple in -y-FelijEj'Eil 
or in frustrated antiferromagnetg23E3. In addition, the 
non-coUinear formulation enables studies of finite tem- 
perature properties of magnetic materials. As the dom- 
inant magnetic excitations at low temperatures are spin 
waves which are non-collinear by nature, it is possible to 
determine the magnon spectra aiid|jJjtiniately the Curie 
temperature from first principlesoE§E3n3. Most of the 
previous work has been done for elements or compounds 
with only one magnetic constituent. We study here sys- 
tems with several magnetic atoms and show how the in- 
teraction between different magnetic sublattices can give 
rise to interesting effects. 



2 



The paper is organized as follows. Some general prop- 
erties of spin spirals are discussed in Sec. || followed by 
the description of the computational scheme in Sec. 111. 
We study the total energy and magnetization with spiral 
magnetic orderings and estimate the Curie temperature 
in Sec. IV and finally we conclude in Sec. 0. 



II. GENERAL PROPERTIES OF SPIN SPIRALS 

The magnetic configuration of an incommensurate spin 
spiral shows the magnetic moments of certain atomic 
planes varying in direction. The variation has a well- 
defined period determined by a wave vector q. When 
the magnetic moment is confined to the lattice sites the 
magnetization M varies as 



M(r„) 



cos(g • r„ + (j)n) sm{9n) 
sm{q ■ r„ + (/)„) sin(6'„) 
cos(6'„) 



(1) 



where polar coordinates are used and m„ is the magnetic 
moment of atom n with a phase 0„ at the position r„. 
Here, we consider only planar spirals, that is On = 7r/2 
which also gives the minimum of the total energy. The 
magnetization of Eq. (Q) is not translationally invariant 
but transforms as 



M(r + i?) = D{q- R)M{r) 



(2) 



where i? is a lattice translation and Z? is a rotation around 
the z-axis. A spin spiral with a magnetization in a gen- 
eral point r in space can be defined as a magnetic con- 
figuration which transforms according to Eq. (^) . As the 
spin spiral describes a spatially rotating magnetization, 
it can be correlated with a frozen magnon. 

Because the spin spiral breaks the translational sym- 
metry, the Bloch theorem is no more valid. Computa- 
tionally, one should use large supercells to obtain total 
energies spin spirals. However, one can define generalized 
translations which contaio-translations in real space and 
rotations in spin spacet§E3. These generalized transla- 
tions leave the magnetic structure invariant and lead to 
a generalized Bloch theorem. Therefore the Bloch spinors 
can still be characterized by a fc-vector in the Brillouin 
zone, and can be written as 



Ah- r 



^-iq-r/2 



Uk{r) 



(3) 



The functions Uk{r) and dk{r) are invariant with respect 
to lattice translations having the same role as for normal 
Bloch functions. Due to this generalized Bloch theorem 
the spin spirals can be studied within the chemical unit 
cell and no large supercells are needed. 

Although the chemical unit cell can be used, the pres- 
ence of the spin spiral lowers the symmetry of the system. 
Only the space group operations that leave invariant the 
wave vector of the spiral remain. When considering the 



general spin space groups, i.e. taking the spin rotations 
into account, the space group operations which reverse 
the spiral vector together with a spiiijptation of n around 
the X-axis are symmetry operationsEHI. 

Basically, the spin spiral relates only the magnetiza- 
tions in the different primitive cells. However, the sym- 
metry properties constrain the magnetization which we 
discuss here in the context of the L2i structure. The 
primitive cell of the L2i structure (one fourth of the cu- 
bic cell shown in Fig. |^) contains four atoms, two Ni, 
one Mn and one Ga or Al atom. In the full cubic sym- 
metry the two Ni atoms are equivalent but this equiv- 
alence can be broken when the spin spiral lowers the 
symmetry of the system. If the spiral wave vector is 
in the [111] direction the two Ni atoms are no longer 
equivalent under space group operations. Considering 
also the spin rotations, the phases </>„ of the two Ni 
magnetizations are opposite as the atoms are related by 
space inversion. If the two Ni atoms are treated equiv- 
alent (when allowed by the spiral symmetry) are con- 
strains for the phases of Ni moments even stronger. If 
the magnetic moments of Ni within the primitive cell are 
M{ri) = micos{(pi) — M{r2), the magnetic moment in 
the neighbouring cell at — ri is M{~ri) = micos{—(j)i). 
On the other hand, the Ni atoms at — ri and at r2 are 
connected by a lattice translation, so that according to 
Eq. (^ M{r2) = rnicos(— 01 + q - R) and one has the re- 
lation 01 = —01 + q - R ior the phase. In order to obtain 
the true minimum energy configuration it may be neces- 
sary to treat the Ni atoms as inequivalent (i.e. lower the 
symmetry of the system) so that the above relation for 
the phase do not have to hold. 



III. COMPUTATIONAL METHOD 



The spin spirals discussed in section |h| are stud- 
ied within the density-functional theory. We use the 
full-potential lineaxiaed augmented-plane-wave method 
(FLAPW) methocCj in an implementation whiciijajlows 
non-coUinear magnetism including spin spiralsElla. In 
addition to the full charge density and to the full poten- 
tial, the full magnetization density is used. The magnetic 
moment is allowed to vary both in magnitude and in di- 
rection inside the atomic spheres as well as in the inter- 
stitial regions. The plane wave cut-off for the basis func- 
tions is RKmax = 9, leading to ^ 350 plane waves with 
the mufhn-tin radii 2.25 a.u. Brillouin zone integrations 
arc carried out with the special point method using 800 k- 
points in the full Brillouin zone and a Fermi broadening of 
0.005 Ry. For the exchange-correlation potential we use 
both the local-spin-density approximation (LSD[A|)EZl and 
the generalized gradient approximation (GGA)l23 which 
we discuss next in more detail. 
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A. LSDA vs GGA 

It has been pointed out that thCjjiat! of GGA is ben- 
eficial in the context of Ni2MnGaO'E£l. Because there 
has been some discussion about the different exchange- 
correlation potentials in the context of non-coUinear mag- 
netism, we present some comparison also here. 

Although there is no global spin quantization axis, one 
can consider at every point of space a local coordinate 
system such that the magnetization at that point is in 
the z-direction. As the LSDA depends only on the mag- 
nitude of the magnetization, the exchange-correlation po- 
tential can be calculated at every point in the local co- 
ordinate system as in the usual coUinear case. The non- 
collinear potential is obtained by rotating back to the 
global frame of reference. On the other hand, the GGA 
depends also on the gradients of magnetization. Because 
the magnetization direction may vary, only projections of 
the magnetization on the local quantization axis are used 
in the standard GGA when evaluating the gradients. If 
the magnetization direction varies slowly this should not 
bring any problems. Some previous work has suggested 
that the disagreements between theory ani. experiment 
are due to projection errors in some casescj. However, 
later work has corroborated that the main issue is not the 
exchange-correlation functional but the actual computa- 
tional method, pointing to the in™o*taHce of all-electron 
and full-magnetization treatmentEl'ESEa. 

We have done all the calculations in this work both 
with LSDA and with GGA. The total energy as a func- 
tion of the spiral wave vector length in Ni2MnGa is shown 
in Fig. for a single direction. One can see that for small 
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FIG. 2: Total energy as a function of the spiral vector q in 
units of 27r/a. Q LSDA, □ GGA. 



q both approximations give similar results. With larger 
q the results differ slightly but the same qualitative be- 



haviour is seen. For the other results presented in the 
following sections the qualitative behaviour is also the 
same for LSDA and GGA, and the quantitative differ- 
ences between the two approximations are even smaller. 
Therefore, only the GGA results are discussed in the fol- 
lowing. 



IV. RESULTS AND DISCUSSION 
A. Total energies 

First, we have studied the possibility of non-coUinear 
ordering by studying the energetics of spiral configura- 
tions. This study also provides information about finite- 
temperature properties. The total energy is calculated 
as a function of the spiral wave vector q, and the wave 
vector is varied along the high symmetry directions [001], 
[110] and [111], q is given in units of l-nja where^fl is the 
theoretical lattice constant of the L2i structureEj. The 
corresponding total energies are shown in Figs. |^ and ^. 

Fig. ^ shows that the variation of total energy in [001] 
and [111] directions is similar in Ni2MnGa and Ni2MnAl 
for all values of q. The lowest energy in all cases is at 
q — Q which is the normal coUinear ferromagnetic config- 
uration. Both materials have small minimum at the an- 
tiferromagnetic configuration at q = (0 1), but at other 
antiferromagnetic configuration oX, q— (0.5 0.5 0.5) there 
are no minima. 




[OOq] [qqq] 

FIG. 3: Total energy as a function of the spiral vector q. Q 
Ni2MnAl, □ Ni2MnGa. 

The energy in the [110] direction is also similar for both 
materials as seen in Fig. ^ Here, the effect of symme- 
try constraints can be seen clearly. If the two Ni atoms 
are equivalent the energy is higher especially around 
q = (0.5 0.5 0). When the magnetic moments of the two 
Ni are allowed to relax independently the energy lowers 
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and the dispersion becomes flat after q = (0.5 0.5 0). 
Near the Brillouin zone boundary at q — (0.75 0.75 0) 
both materials show small energy minima corresponding 
to incommensurate spiral order. At the antiferromag- 
netic configurations at q = (1 1 0) there are no clear 
energy minima even though in the case of Ni2MnAl the 
dispersion is very flat. 
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FIG. 4: Total energy as a function of the spiral vector q. Q 
Ni2MnAl, □ Ni2MnGa. a) Ni atoms are equivalent, b) Ni 
atoms are inequivalent. Vertical line denotes the Brillouin 
zone boundary. 

Generally, the spin spirals are related to magnons 
which allows the estimation of magnon-related proper- 
ties, such as spin stiffness and Curie temperature, from 
the total energies calculated above. The total energy of 
t l^^ nar spin spiral is related to the magnon energy Uq 



4:tiB 
M 



(4) 



where M is the magnetic moment per unit cell. In the 
low q limit the magnon dispersion is quadratic, and one 
defines the spin stiffness constant D as 



(5) 



From the calculated total energies in Figs. I and I we can 
estimate the same spin stiffness for both materials which 
is D = 77 mRy a.u.^. This in good agreement with the 
expecipiental value 79 mRy a.u.^ measured in Ni-Mn-Ga 
filmsN. 

The Curie temperature can be estimated on the basis 
of the Heisenberg model. By mapping the first-principles 
results to the Heisenberg model, the Curie temper|atu*e 
Tq in the random phase approximation is given byEZrEH 



1 



QfiB V 

ksTc M (27r)3 



1 



0J„ 



(6) 



where V is the unit cell volume, and the integration is 
over the Brillouin zone. An estimation can be obtained 
using the quadratic dispersion, Eq. and carrying out 
the integration over a sphere having the same volume as 
the Brillouin zone. This results in 



B^C 



Wqg 

Mtt^D ' 



(7) 



where qd = {6tt^/VY^^. By using the calculated spin 
stiffness constant we obtain Tc=830 K which is clearly 
an overestimate. As seen in Figs. |3| and ^ the disper- 
sion curve E{q) deviates strongly from the quadratic be- 
haviour with larger q. A better estimate can be obtained 
by considering the dispersion quadratic up to some ra- 
dius and constant thereafter. Based on the calculated 
energies in Figs. ^ and ^ the constant is chosen to be 5 
mRy when q > 0.7qd- The Curie temperature obtained 
in this way is Tc—^S5 K which compares well with the 
experimental one 380 K. 



B. Magnetic moments 

In order to obtain a deeper understanding of the energy 
dispersion we next look into the behaviour of magnetiza- 
tion. The magnetic moments averaged over the atomic 
spheres for different q are shown in Figs. || and ^. The 
atomic magnetizations show that within the Mn spheres 
the magnetization is nearly constant and the variation in 
the total magnetization is mainly due to Ni. Also, the 
symmetry consideration of the equivalence of Ni atoms 
has no effect on the Mn moment. This points to a more 
localized character of the magnetic moment of Mn, com- 
pared to a more itinerant character of Ni. Because most 
of the total magnetic moment comes from Mn, these al- 
loys can be considered as localized-moment systems coBr 
sistent with the traditional view for similar materialsES. 
However, despite the relative smallness of its magnetic 
moment, Ni has significant effect for the energetics as dis- 
cussed later on. The differences between Ni2MnGa and 
Ni2MnAl are small: the magnetic moment in Ni2]VLoGa 
is slightly larger as shown already in previous workllj. 

As the magnetic moment in Ni shows a larger vari- 
ation, the behaviour of the Ni moment is analyzed in 
more detail for several directions. The magnetization de- 
creases monotonously both in the [001] and in the [111] 
directions. Differences are at the antiferromagnetic con- 
figurations as the magnetic moment of Ni remains finite 
at q = (1 1 1) but vanishes at q = (0 1). In the 
[110] direction, the behaviour of the Ni moment depends 
strongly on the symmetry as seen in Fig. ^ When the two 
magnetic moments are forced to be the same, the magne- 
tization starts to decrease with increasing q and vanishes 
to zero value at q = (0.5 0.5 0). For larger q values the 
moment shows a small peak before decreasing again to 
zero in the antiferromagnetic state at q = (1 1 0). In the 
case of Ni atoms being inequivalent only a monotonous 
decrease similar to the [001] direction is seen. 
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Because most of the variation in total magnetization 
is due to Ni, it should have a larger effect also on the en- 
ergy dispersion. The importance of Ni can be seen most 
clearly in the [110] direction for the cases of different 
symmetry. The symmetry affects only Ni as seen in the 
behaviour of the magnetization, Fig. As the energy dis- 
persion depends on the symmetry, Fig. ^, the importance 
of Ni is clear. It can be noted also that the energy low- 
ers when the Ni moment increases. Based on the above 
reasoning, Ni should have an effect on the Curie tem- 
perature, which indeed is seen in experiments where the 
increase in Ni content decreases the Curie temperaturety. 
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FIG. 5: Magnetic moments within the atomic spheres as a 
function of the spiral vector q. Q Ni2MnAl, □ Ni2MnGa. a) 
Ni atoms are equivalent, b) Ni atoms are inequivalent 
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FIG. 6: Magnetic moments within the Ni sphere as a function 
of the spiral vector q. Q Ni2MnAl, □ Ni2MnGa. 



The variation of the Ni moment can be understood by 



considering symmetry arguments and the coordination 
around Ni atoms. In the [001] direction two of the four 
Mn atoms neighbouring Ni have the same magnetization 
direction in the spiral and the other two have different 
direction, as shown schematically in Fig. |7|(a). The mag- 
netization in Ni favours ferromagnetic alignment with the 
neighbouring Mn moments so that part of the Ni moment 
can be thought to align with one group of the Mn neigh- 
bours and part with the other group. The total moment 
within the atomic sphere is average of these two parts 
and the Ni moment decreases when the angle between the 
Mn moments increases. In the antiferromagnetic config- 
uration there is a complete frustration of the Ni atoms 
which results in the zero average magnetization within 
the sphere. For the [111] direction shown in Fig. 0(b) 
one group contains three Mn atoms and the other group 
only one. Therefore the variation of the average moment 
in the Ni sphere is smaller and the moment remains finite 
in the antiferromagnetic configuration. 





FIG. 7: Schematic view of the magnetic moments in the 
nearest neighbour Mn atoms of Ni at a) q = (0 1) b) 
q = (0.5 0.5 0.5). 

In the [110] direction the situation is more complex es- 
pecially when the two Ni atoms are treated as equivalent. 
In the antiferromagnetic configuration the coordination 
is similar to the case of the [001] direction. There are two 
groups of neighbouring Mn with antiparallel magnetiza- 
tion, and the frustration leads to zero average moment 
within the Ni sphere. The Ni moment is, however, zero 
also a.t q — (0.5 0.5 0). At this point there are three 
groups of equivalent Mn neighbours. One group con- 
tains two Mn atoms and the other groups contain one 
Mn atom. The magnetic moments of single Mn atoms 
are antiparallel to each other and have 90° angle with re- 
spect to the moments in the group of the two Mn atoms. 
The other equivalent Ni atom has three similar groups 
of neighbouring Mn as the first Ni. The important point 
is that the moments in the group with two Mn atoms 
are antiparallel to those in the corresponding group of 
the first Ni, as seen in Fig. ||. There is now frustra- 
tion for Ni, but only when both equivalent Ni atoms and 
their neighbours are taken into account. This frustration 
causes the magnetic moment around Ni to vanish com- 
pletely in contrast to the antiferromagnetic case, where 
a small moment remains near Ni but averages to zero. 
When the Ni are inequivalent, they can relax according 
the local environment so that a finite moment can remain 
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&tq = (0.5 0.5 0) 






FIG. 8: Schematic view of the magnetic moments in the 
nearest neighbour Mn atoms of the two equivalent Ni at 
q = (0.5 0.5 0). 

An example of the magnetization density for the case 
where finite magnetic moments near the Ni atom average 
to zero is seen in Fig. ^. Here the magnetization direction 
can change its sign within the atomic sphere. This finding 
shows the importance of the full magnetization treatment 
when dealing with several magnetic sublattices. 



FIG. 9: Magnetization density around Ni in the (001) plane 
with q = (1 1 0). The width and the height of the area are 
2.5 a.u. while magnetization is in arbitrary units. 



V. CONCLUSIONS 

We have studied non-coUinear magnetic configurations 
in the ternary alloys Ni2MnGa and Ni2MnAl with first- 
principles calculations. The calculations show that the 
magnetic properties are similar for both materials. The 
ferromagnetic configuration is the ground state in the L2i 



structure, so that the experimentally observed antiferro- 
magnetism of Ni2MnAl is related to structural disorder. 
The calculated total energies are used to estimate the 
spin stiffness constant and the Curie temperature which 
are in good agreement with the experiments. The simi- 
larity in the energy dispersion for both materials suggests 
that the Curie temperatures should be also similar. In 
the [110] direction Ni2MnAl has higher energy, so that 
the Curie temperature should be a slightly higher. 

The variation of the magnetic moment in the spirals 
shows that the Mn moment is nearly constant while the 
Ni moment varies strongly. The symmetry of the spin 
spiral constrains the direction of magnetization, and as 
Ni favours ferromagnetic coupling with Mn, there can be 
frustration at certain wave vectors resulting in the van- 
ishing of the magnetic moment near the Ni sites. It is also 
shown how there can be strong variation in the direction 
of the magnetization near the atomic sites which points 
to the relevance of the full magnetization treatment. 

Some conclusions can be made concerning the role of 
the constituent atoms for the magnetic properties. As 
the magnetic moment of Ni varies strongly and its sym- 
metry affects the energy considerably, Ni has probably 
a strong effect on the energy dispersion especially when 
larger wave vectors are involved. Therefore Ni also in- 
fluences the Curie temperature. If one assumes that the 
spin stiffness is mainly due to Mn, and the lowering of 
the energy with larger wave vectors due to Ni, Ni low- 
ers the Curie temperature from 830 K to 485 K within 
the present approximations. As the increase in the Ni 
moment decreases the energy it is suggested that in or- 
der to increase the Curie temperature one should replace 
some Ni, perhaps a little counter-intuitively, with some 
non-magnetic element, for example Cu. Further experi- 
ments should clarify these issues and confirm the above 
suggestions. 
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